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NAIE> Preface 
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¢ This tutorial introduces terminology and concepts 
used in the later SPICE tutorials. 


¢ Some of this material is more difficult than what 
follows in later presentations. 


— Acomplete understanding of this material is not essential in 
order to use SPICE. 


¢ Still, we think this information may be helpful, 
SO... on we go! 
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Time 
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¢ An epoch is an instant in time specified by some 
singular event 
— Passage of a star across your zenith meridian 
— Eclipse of a spacecraft signal as it passes behind a solid body 


¢ Clocks 


— Clocks count epochs specified by events such as: “regular” 
oscillations of a pendulum, quartz crystal, or electromagnetic 
radiation from a specified source, measured from an agreed 
upon reference epoch. 


— Careful specification of epochs using clocks requires reference 
to the particular clock and the location of that clock. 


¢ Time Systems 


— Are agreed upon standards for “naming” epochs, measuring 
time, and synchronizing clocks 
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Atomic Time and UTC 
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¢ International Atomic Time (TAI) 
— Statistical time scale 


» Based on data from ~200 atomic clocks in over 50 national 
laboratories 


— Maintained by Bureau International des Poids et Mesures (BIPM) 
— Unit is the SI (System International) second 


» duration of 9192631770 periods of the radiation corresponding to 
the transition between two hyperfine levels of the ground state of 
the cesium 133 atom 


— TAI is expressed as a count of atomic seconds past the astronomically 
determined instant of midnight 1 Jan 1958 00:00:00 


¢ Coordinated Universal Time (UTC) 
— Civil Time at Greenwich England (~GMT) 
— Usual Calendar Formats plus Hour:Minute:Second.fraction 


¢ Relationship between TAI and UTC 
— UTC + 10 seconds + number of leap seconds = TAI 
» Valid only after Jan 01, 1972 
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Astronomical Time — UT1 
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¢ Astronomical Time (UT1) is an hour 
representation of the angle between the 
Greenwich zenith meridian and the location of the 
“computed mean sun.” 


¢ Used prior to atomic time for civil time keeping 
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Tying UTC to Earth’s Rotation (UT1) 
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Ideally, UTC noon and astronomical noon at Greenwich (UT1) should occur 
simultaneously. However, the earth’s rotation is not uniform. Eventually, 
UTC noon and astronomical noon at Greenwich get out of sync. 


UTC Ka 


(ro 


Tye 3 . Pg 


When the mismatch becomes greater than 0.9°atomic 
seconds, a “leap second” is added to (or removed from) 
the end of a designated UTC day—normally either June 30 
or December 31. 


The variations in the earth’s rotation that cause leap 
seconds to be needed are not predictable. 
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Leapseconds (+ and -) 
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- “Normal” sequence of ° Sequence with a 


UTC time tags Positive Leapsecond 
-— 1998 Dec 31 23:59:58.0 - 1998 Dec 31 23:59:58.0 
-— 1998 Dec 31 23:59:59.0 - 1998 Dec 31 23:59:59.0 
-— 1999 Jan 01 00:00:00.0 - 1998 Dec 31 23:59:60.0 
-— 1999 Jan 01 00:00:01.0 - 1999 Jan 01 00:00:00.0 
-— 1999 Jan 01 00:00:01.0 
Leap seconds complicate the task of ‘ Sequence with a 
finding the duration between two ‘ 
UTC epochs: Negative Leapsecond 
-— 1998 Dec 31 23:59:57.0 
¢ You need to know when past leap ~ 1998 Dec 31 23:59:58.0 
seconds occurred to compute 
durations defined by pairs of past - 1999 Jan 01 00:00:00.0 
UTC epochs. —~ 1999 Jan 01 00:00:01.0 


¢ Durations defined by pairs of future 
UTC epochs are indeterminate if leap 
seconds could occur in the interim. 
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Barycentric Dynamical Time 
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¢ Barycentric Dynamical Time (TDB) and Ephemeris 
Time (ET) are synonyms in SPICE documentation. 


¢ TDB is 
— a mathematical ideal used in the equations of motion. 


— used as the independent time variable for many SPICE 
subroutine interfaces. 


— related to Barycentric Coordinate Time (TCB) by an offset and 
a scale factor. 


— TDB advances on average at very close to the same rate as 
TAI---the difference is nearly periodic. 
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NAIE> Terrestrial Dynamical Time 
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¢ Terrestrial Dynamical Time (TDT) 


— TDT is the Ideal Time (proper time) on Earth at sea level 
— TDT = TAI + 32.184 seconds 


— The IAU has adopted the name “Terrestrial Time” (TT) 
» But this is called TDT throughout SPICE documentation 
¢ TDB and TDT have nearly the same reference 
epoch (approximately 1 Jan 2000, 12:00:00 at 
Greenwich England), called “J2000.” 
¢ TDB and TDT advance at different rates. 
— Variations are small: ~ 1.6 milliseconds 


— Variations are almost periodic with a period of 1 sidereal year 
(to first order) 


— Variations are due to relativistic effects 
» TDB = TDT + 0.001657 sin( E + 0.01671sin(E) ) 
¢ Use of TDT in the SPICE system is quite limited. 
— SCLK kernels 
— Duration computations involving UTC 


Fundamental Concepts 


Offsets between Time Systems 
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Difference between seconds past J2000 in a given time system 
and TDB seconds past J2000 TDB. Systems used for 
comparison are TDT, TAI, and UTC: 


TDT-TDB, TAI-TDB, UTC-TDB vs TDB 
10 T T T T T T 


|| _ TDT-TDB_ 


Magnified by 1000 


Actual amplitude of 
periodic term = = SEG SREEERESSSRECSECVEREEAB ECHO MaldaleaisibnSaicibipindesinn TN 
~1.6e-3 seconds 


ET DB. scwisscsrsssssnstannsoosrunsivare ntnaeiecigll kecetcousetester eusssaoamnai | 
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Fundamental Concepts 11 


NAIE> Spacecraft Clocks 
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¢ Spacecraft have onboard clocks to control 
scheduling of observations, maneuvers, attitude 
adjustments, etc. 


¢ Used to time stamp data 
¢ Fundamental unit of time is the “tick” 


— Smallest increment possible for a spacecraft clock 
— Nominal tick duration is spacecraft clock dependent 


¢ Spacecraft clock time is a count of ticks since 
some reference tick. 


¢ The duration of the tick drifts with respect to other 
time systems because spacecraft clocks are not 
very stable 
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More about Spacecraft Clocks 
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¢ SCLK string formats vary from one spacecraft clock to the 


next. 

— Cassini: Maximum reading for partition 1 = 1/4294967295.255 
» Partition number: 1 
» Seconds: 4294967295 
» Ticks (for Cassini, unit = 1/256 second): 255 

— Galileo: Maximum reading for partition 1 = 1/16777215:90:09:07 
» Partition number: 1 
» "RIM" count (unit = 60 2/3 seconds): 16777215 
» "Mod 91" count (unit = 2/3 second): 90 
» "RTI" count (unit = 1/15 second): 9 
» "Mod 8" count (unit = 1/120 second): 7 


¢ The format of spacecraft clock and the relationship between 
tick count and other time systems (usually UTC) is captured 
ina SPICE SCLK kernel 
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SPICE Definitions: 
Reference Frames & Coordinate Systems 
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¢ A reference frame is an ordered set of three 
mutually orthogonal (possibly time dependent) unit- 
length direction vectors, coupled with a location 
called the frame’s “center” or “origin.” 


— SPICE documentation frequently uses the shorthand 
“frame” instead of “reference frame.” 


— The ordered set of axes of a reference frame is also called a 
“basis.” 
¢ A coordinate system specifies the method of 
locating a point within a reference frame. 


iZ ; 


Two examples of coordinate systems 


-———--s 


Cartesian coordinates Spherical coordinates 
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Reference Frame Center 
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¢ A reference frame’s center is an ephemeris object whose 
location is coincident with the origin (0, 0, 0) of the frame. 
— The center of the IAU_<body> frame is center of mass of <body>. 
— The center of any inertial frame is (in SPICE) the solar system barycenter. 
» True even for frames naturally associated with accelerated bodies, 
such as MARSIAU. 
¢ A frame’s center plays little role in specification of states 
— Origin cancels out when doing vector arithmetic 


» Whether positions of objects A and B are specified relative to centers 
C1 or C2 makes no difference: 


(A-C1)-(B-C1)=(A-C2)-(B-C2)=A-B 


— But the center *is* used in computing light time to centers of non-inertial 
frames 


» When the aberration-corrected state of Titan as seen from the Cassini 
orbiter is computed in the body-fixed IAU_ Titan frame, light time is 
computed from Titan’s center to the Cassini orbiter, and this light 
time is used to correct both the state and orientation of Titan. 


» (Light time and aberration corrections are discussed later on.) 
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Types of Reference Frames 


Navigation and Ancillary Information Facility 


¢ Inertial 
— Non-rotating 
» With respect to fixed stars 
— Non-accelerating origin 
» Velocity is typically non-zero; acceleration is negligible 
— Examples: 
» J2000 (also called ICRF), B1950 


¢ Non-Inertial 
— Examples 

» Body-fixed 
¢ Centered at body center 
¢ Topocentric 

» Instrument 

» Dynamic frames 
¢ For example, frames defined by time-dependent vectors 
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J2000 Frame 
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¢ The J2000* (aka EME2000) frame definition is based 
on the earth’s equator and equinox, determined 
from observations of planetary motions, plus other 
data. 


Z is normal to the mean equator of date at epoch J2000 
Zj2000 | TDB, which is approximately Earth’s spin axis orientation 

at that epoch. (J2000 TDB is 

2000 JAN 01 12:00:00 TDB, or JD 2451545.0 TDB). 


Equatorial plane 

Plane normal to the earth’s 

apiiexier <— Ecliptic plane 
Plane defined by movement 
of the earth around the sun 


a 
\ a Y 42000 
ON Y =Z cross X 


Intersection of 
X 52000 equatorial and 
ecliptic planes, 
called vernal 
equinox 


*Caution: The name “J2000” is also used to refer to the zero epoch of the ephemeris time system (ET, also known as TDB). 
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The ICRF Frame 
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¢ The ICRF* frame is defined by the adopted 
locations of 295 extragalactic radio sources 


Solar System Barycenter 


*ICRF = International Celestial Reference Frame 


‘ ‘ j : R 
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J2000 vs. ICRF 
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¢ The realization of ICRF was made to coincide almost 
exactly with the J2000 frame. 
— The difference is very small—a rotation of less than 0.1 arc second 
— These two frame names are treated synonymously in SPICE 


» In reality, any SPICE data said to be referenced to the J2000 
frame are actually referenced to the ICRF frame 


» For historical and backwards compatibility reasons, only the 
name “J2000” is recognized by SPICE software as a frame 
name-not “ICRF” 


¢ No transformation is required to convert SPICE state 


vectors or orientation data from the J2000 frame to 
the ICRF. 


ICRF = International Celestial Reference Frame, defined by the IERS 
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Body-fixed Frames 


Navigation and Ancillary Information Facility 


¢ Body-fixed frames are tied to a named 
body and rotate with it 


— The most common names, those for planets, 
satellites, the sun, and a few asteroids and comets, 
are hard-coded in SPICE software 


» Frame name style is “IAU_body name” 
+ Examples: IAU_MARS, IAU_SATURN 
» To see all such names, see: 
« Frames Required Reading document, or 
« Latest generic PCK file 
— The rotation state (the orientation at time T) is 
usually determined using a SPICE text PCK 
containing data published by the IAU xX 
— The earth and moon are special cases! 


» IAU_EARTH and IAU_MOON both exist but 
generally should NOT be used 


» See the SPICE tutorial named “lunar- 
earth_pck-fk” for the best frames to be used 
for those bodies 


— On very rare occasions a CK is used to provide a 
body’s rotation state 
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Body-fixed 
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Topocentric Frames 
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¢ Topocentric frames are Z points “up” 
attached to a surface Position Vector 


¢ Z-axis is parallel to the 
gravity gradient or 
orthogonal to reference 
spheroid 


x points North 


Azimuth (increases 
in clockwise 
direction, measured 
ee ae from +x axis) 


Elevation (angle between pdeodanal 
vector and x-y plane) reiecienck 
position vector 
One example of a topocentric frame. There onto x-y plane 
are other types of topocentric frames: for 
example, the z-axis could point down, the x- 
axis North, and the y-axis East. 


y points West 
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Coordinate System Conventions - 1 
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¢ Planetocentric coordinate systems (also known as latitudinal) 
— For planets and their satellites the +Z axis (+90 


LAT) always points to the north side of the 
invariable plane (the plane whose normal 
vector is the angular momentum vector of the ‘ 
solar system) 
» Planetocentric longitude increases positively 
eastward 


» Planetocentric latitude increases positively 
northward 


Dwarf planets*, asteroids and comets spin in 
the right hand sense about their “positive 


pole.’ as: | 
» What the IAU now calls the “positive pole” is Y 


still referred to as the “north pole” in SPICE 
documentation. 


» The “positive pole” may point above or below 
the invariable plane of the solar system (see 
above). 


» This revision by the [AU Working Group (2006) 
inverts what had been the direction of the north 
pole for Pluto, Charon and Ida. 

Planetocentric routines are RECLAT/ LATREC, 
RADREC/RECRAD, DRDLAT/DLATDR 


xX 


*The dwarf planets are: Ceres, Eris, Haumea, Makemake, Pluto 
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Coordinate System Conventions - 2 
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¢ Planetodetic coordinate systems 


— Planetodetic longitude is 
planetocentric longitude 


» Increases positively eastward 
— Planetodetic latitude 
» Tied to a reference ellipsoid 


» Fora point on a reference 
ellipsoid, the angle measured 
from the X-Y plane to the 


surface normal at the point of Vass J 
interest. For other points, equals es — : 
latitude at the nearest point on a, 
the reference ellipsoid 
» Latitude increases positively in 
the +Z direction 
¢ +Z is defined the same as for 
planetocentric coordinates 
— Planetodetic routines are 
RECGEO/GEOREC, 
DRDGEO/DGEODR 
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Coordinate System Conventions - 3 
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¢ Planetographic coordinate systems 


— For planet and satellite planetographic 
coordinate systems: 


» Planetographic latitude is planetodetic z 
latitude Spin 


direction 


» Planetographic longitude is usually 
defined such that the sub-observer 
longitude increases with time as seen by 
a distant, fixed observer in an inertial 
reference frame 


¢ The earth, moon and sun are exceptions; 
planetographic longitude is positive east by 
efau 


* Planetographic routines are PGRREC/RECPGR 


— For dwarf planets, asteroids and comets: x 


» There are Was inconsistent 
standards! (USNO, IAU, PDS) 


» NAIF strongly suggests you use only | 
planetocentric or planetodetic 
eoote nae 


Planetocentric routines are RECLAT/ LATREC, 
RADREC/RECRAD, DRDLAT/DLATDR 
- Planetodetic routines are RECGEO/GEOREC, Distant Observer 
Fundamental Concepts DRDGEO/DGEODR 
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State Vectors 
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¢ The state of an object is its position and velocity 
relative to a second object 


— In SPICE, these objects are often referred to as “target” and 
“observer” or “center” 


— E.g. Saturn relative to Saturn barycenter; Titan relative to Huygens 
probe 


— A state is always given in the TDB time system (also known as ET) 
¢ In the SPK subsystem a state is a six dimensional 
vector 
— First three components are Cartesian position: x, y, z 


— Second three components are Cartesian velocity: dx/dt, dy/dt, dz/dt 
— Units are km, km/sec 


¢ A state is specified relative to a reference frame 
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CNAIE> Transforming States 
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¢ To perform algebraic operations on states they 
must be in the same reference frame. 


¢ Position-only frame transformations require only 
a rotation® matrix given as a function of time. 
» Pp (f) = Ra to p(t) Pa(t) 
¢ Position and velocity frame transformations 
require that we differentiate the above equation 
» AP (f) dt = ARatoa(t/dt Pa(t) + Ratoa(t) d Pa(é/dt 
¢ We can use a 6x6 matrix to combine these two 
transformations into a single equation 
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CNAIE> Transforming States 
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S_a(t) = Tato a(t)Sa(t) 


where 


(" ) 
Si(f) = 
AP;(t)/dt i=AorB 


and 


Ra to p(t) 0 
Tato p(t) = ( ) 
dRa to p(t)/dt Ra to p(t) 


The SPICELIB routines SXFORM and PXFORM return state 
transformation and position transformation matrices respectively. 
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Aberration Corrections: Introduction 
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¢ Within the SPICE system, “aberration corrections” 
are adjustments made to state vectors and time- 
dependent reference frames to accurately reflect 
the apparent—as opposed to the actual-—state and 
attitude of a target object as seen from a specified 
observer at a specified time. 


— Actual, uncorrected states from an ephemeris are called 
“geometric” states. 


— When computing state vectors, SPICE users may request 
geometric or aberration-corrected states. 
¢ Aberration corrections are needed to accurately 
answer questions such as: 


— In which direction must a remote sensing instrument be pointed 
to observe a target of interest? 


— For a given pointing direction and observation time, what target 
body surface location would be observed by a remote sensing 
instrument? 


— In which direction must an antenna be pointed to transmit a 
signal to a specified target? 
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Computing Aberration-corrected States 
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¢ In order to compute an aberration-corrected state, 
the state vectors used in the computation must all 
be chained back to the Solar System Barycenter 


Saturn 


Cassini spacecraft 
during interplanetary 
cruise 


The aberration-corrected red ea System 
vector cannot be computed Barycenter 


unless all three green ‘ + ‘ 
vectors are available. . (*) 
He 


Sun 
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Example: Predicted vs Actual Photo 
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We compare the predicted appearance of a photograph from an optical camera 
against the actual photograph. We show three predictions derived using different 
aberration corrections: NONE, LT ("light time only"), and LT+S (“light time plus 
stellar aberration"). 


For each prediction, we use red overlays to indicate the expected location in the 
photo of the images of an extended target body (for example, a natural satellite), of 
features on the surface of the target body, and of a star. 


@ lee 


NONE LT+S 
Predicted images using Predicted images using light Predicted images using light 
uncorrected target position time-corrected target time and stellar aberration- 
and orientation and position and orientation and =‘ corrected target position, light 
uncorrected star direction uncorrected star direction time-corrected target 
vector vector orientation, and stellar 


aberration-corrected star 


direction vector 
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Prediction Without Corrections 
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Predicted target body, surface feature, 
and star image locations (in red). Actual 
locations are in blue and black. 


Any Body motion 


Actual == relative to solar 
star system barycenter 


image 


Using geometric target positions, target images in photos or other remote-sensing 
observations don’t appear at their predicted locations. 
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Light Time Corrections 
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Target’s 
velocity Target’s angular 
velocity vector 


Target’s (geometric) position 


and orientation at time ET Target’s position and 


orientation at time ET-LT 
(light time corrected 


Note: The angular separation of ve ‘ 
' position and orientation) 


the geometric and light time 

corrected positions as seen by 
the observer has been 
exaggerated for readability. 


Observer’s position at time ET 


At time ET, the observer’s camera records photons emitted from the target at 
time ET-LT, where LT is the one-way light time. The camera "sees" the 
target's position and orientation at ET-LT. 
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Prediction Using Light Time Corrections 
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Predicted target body, surface feature, 
and star image locations (in red). Actual 
locations are in blue and black. 


sy Body motion 


relative to solar 


Actual system barycenter 


star 
image 


Using the light time corrected target position and orientation, the predicted 
locations in the photo of the target image and surface features have changed, but 
the accuracy of the prediction may still be poor. 
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Stellar Aberration Correction 
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Target’s apparent position 


Target’s position and orientation and orientation 


at time ET. This position is called 
the geometric or uncorrected 
position. Target’s position and 
orientation at time ET-LT 
(light time corrected 


position and orientation) 


Note: The angular separation of 
the geometric, light time 
corrected, and apparent 
positions as seen by the 
observer has been exaggerated 
for readability. 


Velocity component 
orthogonal to 

observer-target 
position 
To the observer, light appears to travel this path 
Observer’s velocit 


y 
relative to Solar System QObserver’s position at time ET 
Barycenter 


At time ET, the observer’s camera records photons emitted from the target at 
time ET-LT, where LT is the one-way light time. 


The vector from the observer at ET to the location of the target at ET-LT is 
displaced by a physical phenomenon called stellar aberration. The displaced 
vector yields the apparent position of the target. 


Fundamental Concepts 


Prediction Using "LT+S" Corrections 
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Actual 
star 
image 


+ 


Actual 
body 
image 


Predicted target body, surface feature, and star 
image locations (in red) 


Using the light time and stellar aberration-corrected target position, light time- 
corrected target orientation, and stellar aberration-corrected star direction, we 
obtain a significantly improved image location predictions. 


Remaining prediction errors may be due to, among other causes, pointing 
error, spacecraft and target ephemeris errors, and timing errors. 
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Newtonian Stellar Aberration Effect 
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ina emission 
4 from Target 


t ET 
Photon path length = c * LT at ET - LT 


“V” = observer velocity component orthogonal to radiation path. 


Observer path 
length = V*LT 


Observer at 
ET-LT 


We ignore parallel velocity component because it has negligible 
effect on photon path geometry. 


Photon 
emission 
from Target 
at ET - LT 


Component of photon path 
orthogonal to V. Length ~= c * LT 


Apparent direction to target 
a ; 
Photon path in observer frame 


Stellar aberration angle dtheta 
sin(dtheta) = v/c +O(v2/c2) 


Component of 
photon path parallel 
to V. Length ~= 
Y, 


-Observer 
velocity 


Camera 


component 
orthogonal to 
radiation path 


Photon 
reception 
at ET 
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NAIE> Effect of Aberration Corrections - 1 
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¢ Angular offsets between corrected and 
uncorrected position vectors over the time span 
2004 Jan 1--2005 Jan1 
— Mars as seen from MEX: 
» LT+S vs NONE: — .0002 to .0008 degrees 
» LT vs NONE: .0006 to .0047 degrees 
— Earth as seen from MEX: 
» LT+S vs NONE: .0035 to .0106 degrees 
» LT vs NONE: .0000 to .0057 degrees 
— MEX as seen from Earth: 
» LT+S vs NONE: _ .0035 to .0104 degrees 
» LT vs NONE: .0033 to .0048 degrees 
— Sun as seen from Mars: 
» LT+S vs NONE: _ .0042 to .0047 degrees 
» LT vs NONE: .0000 to .0000 degrees 
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Effect of Aberration Corrections - 2 
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Angular offsets between corrected and uncorrected position 
vectors over the time span 2004 Jan 1--2008 Jan1 
— Saturn as seen from CASSINI: 
» LT+S vs NONE: .0000 to .0058 degrees 
» LT vs NONE: .0001 to .0019 degrees 
— Titan as seen from CASSINI: 
» LT+S vs NONE: .0000 to .0057 degrees 
» LT vs NONE: .0000 to .0030 degrees 
— Earth as seen from CASSINI: 
» LT+S vs NONE: .0000 to .0107 degrees 
» LT vs NONE: .0000 to .0058 degrees 
— CASSINI as seen from Earth: 
» LT+S vs NONE: .0000 to .0107 degrees 
» LT vs NONE: .0000 to .0059 degrees 
— Sun as seen from CASSINI: 
» LT+S vs NONE: .0000 to .0059 degrees 
» LT vs NONE: .0000 to .0000 degrees 
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